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Abstract 


The  fuel  2-azido-N,N-dimethylethanaiiune  (DMAZ)  has  shown  considerable  promise  as  a 
replacement  for  hydrazine-based  fuels  in  hypergolic  propulsion  systems,  but  the  ignition  delays 
observed  for  DMAZ-inhibited  red  fuming  nitric  acid  (IRFNA)  systems  are  longer  than  those  for 
monomethylhydrazine-IRFNA  systems.  This  report  considers  the  potential  of 
2-azidodocycloaIkanamine-based  fuels  for  addressing  this  issue.  Such  molecules  have  two 
stereochemically  distinct  isomers,  one  of  which  prevents  the  azide  group  from  “shielding”  the 
amine  lone  pair  electrons  from  proton  attack.  The  other  promotes  such  shielding.  Since 
shielding  likely  influences  the  manner  in  which  nitric  acid  reacts  with  amine-azide  fuels,  the 
ignition  delays  for  the  isomers  may  be  different,  and  it  is  possible  that  one  of  the  isomers  will 
yield  shorter  ignition  delays  than  DMAZ.  To  support  consideration  of  the  synthesis  and 
development  of  2-azido-N,N-dimethylcyclopropanamine  (ADMCPA)  equilibrium  structures  are 
identified  and  their  gas-phase  heats  of  formations  estimated  via  density  functional  theory-based 
calculations.  In  addition,  ADMCPA  isomer  condensed-phase  densities  and  heats  of  vaporization 
are  estimated  via  molecular  dynamics  simulations  performed  with  halogen  analogs.  Combined 
and  employed  to  predict  ballistic  properties,  the  computational  results  indicate  that  the 
ADMCPA  isomers  merit  further  investigation  as  hypergols.  A  suggestion  for  synthesizing 
CZ5- ADMCPA  is  also  proffered. 
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1.  Introduction 


Providing  an  extremely  reliable  basis  on  which  to  design  intermittent  and/or  variable  thrust 
propulsion  systems,  hypergolic  liquid  (or  gel)  fuel-oxidizer  combinations — ^i.e.,  those  that  react 
spontaneously  upon  mixing — are  employed  throughout  the  Department  of  Defense  (DOD)  and 
National  Aeronautics  and  Space  Administration  (NASA)  for  aviation  and  missile  applications. 
One  of  the  drawbacks  to  many  “standard”  hypergolic  fuel-oxidizer  combinations,  however,  is 
that  the  fuels  are  derived  from  hydrazine,  monomethylhydrazine  (MMH),  and/or  unsymmetrical 
dimethylhydrazine  (UDMH) — all  of  which  are  acutely  toxic  and  suspected  carcinogens  [1]. 
Searching  for  alternatives  to  hydrazine-based  fuels,  the  Army  (via  the  U.S.  Army  Aviation  and 
Missile  Command  [AMCOM]  and  in  collaboration  with  other  DOD  and  NASA  organizations)  is 
investigating  the  use  of  various  2-azidoethanamines  in  such  applications.  Collectively  referred 
to  as  Competitive  Impulse  Non-Carcinogenic  Hypergols  (CINCH)  [2],  the  formulation  within 
this  class  of  compounds  receiving  the  most  attention  is  2-azido-N,N-dimethylethanamine 
((CH3)2NCH2CH2N3).  Also  referred  to  as  DMAZ,  this  fuel  has  been  foimd  to  yield  performance 
competitive  with  Aerozine  50  (a  50/50  mixture  of  hydrazine  and  UDMH)  in  inhibited  red  fuming 
nitric  acid  (IRFNA)  oxidized  experimental  systems.  DMAZ-IRFNA  systems  do  not,  however, 
meet  “ignition  delay’  standards  set  by  MMH-IRFNA  systems,  and  this  shortcoming  may 
negatively  impact  the  Army’s  ability  to  field  this  alternative. 

Based  on  traditional  chemical  knowledge  and  testing  performed  for  AMCOM,  it  was  concluded 
that  a  rate-limiting  step  in  DMAZ-IRFNA  ignition  might  be  associated  with  an  initial  chemical 
reaction — namely,  proton  transfer  from  nitric  acid  to  DMAZ’s  amine  nitrogen  [3].  As  a  result, 
one  of  the  approaches  to  addressing  the  ignition  delay  issue  has  been  to  synthesize  differently 
substituted  2-azidoethanamines  with  the  hope  of  enhancing  the  amine’s  reactivity/basicity.  Two 
compounds  synthesized  subsequent  to  the  identification  of  DMAZ  as  a  promising  hypergol  are 
2-azido-N-methylethanamine  (MMAZ)  and  2-azido-N-cyclopropylethanamine  (CPAZ),  and  a 
comparison  of  their  attributes  is  instructive.  MMAZ,  which  has  a  hydrogen  atom  in  lieu  of  one 
of  the  methyl  groups  in  DMAZ,  is  not  hypergolic,  while  CPAZ,  which  has  a  cyclopropyl  group 
in  lieu  of  the  methyl  group  in  MMAZ,  is.  However,  the  hypothesis  that  the  ignition  delays  for 
these  fuels  are,  in  fact,  due  to  differences  in  the  basicity  of  their  amine  nitrogen  remains  to  be 
proven. 

In  support  of  the  experimentally  based  research  being  pursued  through  AMCOM,  the  U.S.  Army 
Research  Laboratory  (ARL)  is  employing  computational  quantum  chemistry  in  an  attempt  to 
identify  the  differences  between  fuels  that  correlate  with  ignition  delay,  the  initial  focus  of  this 
effort  being  the  characterization  of  fuel-nitric  acid  proton  transfer  reactions.  Toward 
characterizing  DMAZ-nitric  acid  reactions,  a  study  to  establish  the  nature  of  DMAZ  equilibrium 
stmctures  was  conducted  [4],  with  12  conformers  being  identified  via  density  functional  theory 
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(DFT).  One  of  the  interesting  findings  of  that  study  was  that  in  DMAZ’s  lowest  energy  structure 
(Figure  1),  the  azido  group  will  shield  the  amine  lone  pair  electrons  from  proton  attack.  Given 
the  importance  of  proton  transfer  reactions  inferred  from  experiments,  this  finding  raised  the 
possibility  of  molecular  geometry  having  an  influence  in  the  ignition  process. 


Figure  1 .  The  lowest  energy  DMAZ  structure  and  a  mapping  of  the 

electrostatic  potential  on  its  0.004  atomic  unit  (a.u.)  electron 
isodensity  sinface.  Dark  blue  regions  of  the  potential 
correspond  to  nucleo|diilic  (lone  pair  electron)  sites. 

While  considering  means  to  explore  this  possibility,  doubts  about  the  value  of  a  computational 
study  arose,  the  concern  being  that  even  if  a  deleterious  relationship  between  shielding  and 
ignition  delay  could  be  established,  an  inability  to  prevent  shielding  would  make  the  knowledge 
moot.  This  prompted  the  author  to  consider  alternatives  to  2-azidoethanamines.  Informed  that 
the  introduction  of  a  cyclic  structure  can  “lock  in”  a  preferred  orientation  between  functional 
groups,  the  author  sketched  out  a  number  of  notional  molecules  based  on  the  approach.  Of  the 
molecules  sketched,  2-azidocycloalkanamines,  the  trans  form  of  which  prevents  the  azide  group 
from  coming  in  close  proximity  to  the  amine  lone  pair,  were  intuitively  appealing. 
Stoichiometrically,  2-azido-N,N-dimethylcyclopropanamine  (ADMCPA),  2-azido-N- 
methylcyclobutanamine  (AMCBA),  and  2-azidocyclopentanamine  (ACPA)  differ  from  DMAZ 
by  only  one  carbon  atom,  (C5H10N4  vs.  C4H10N4),  and,  lUce  DMAZ,  the  amine  and  azide  groups 
are  separated  by  two  carbon  atoms.  Moreover,  like  CPAZ,  which  also  has  C5H10N4 
stoichiometry,  the  cycloalkyl  group  should  increase  the  reactivity  of  the  amine  lone  pair. 
Therefore,  if  shielding  of  these  electrons  by  the  azido  group  is  indeed  a  rate-determining 
mechanism  in  the  ignition  process,  the  cis  and  trans  forms  of  a  given  azidocycloalkanamine  will 
have  different  ignition  delays,  and  it  is  possible  that  one  of  the  forms  will  yield  shorter  ignition 
delays  than  DMAZ  or  CPAZ.  AMCBA  and  ACPA  will  presumably  be  less  reactive  than 
ADMCPA,  but  might  be  worth  synthesizing  if  ADMCPA  proves  to  react  too  violently — a. 
problem  encountered  in  CINCH  molecules  synthesized  with  more  than  one  azide  functional 
group  [2, 3]. 


2 


It  is,  of  course,  a  nontrivial  exercise  to  experimentally  synthesize  and  characterize  a  new 
compound,  and  the  intuition  of  a  computational  chemist  is  hardly  adequate  justification  for 
imdertaking  such  a  project.  Therefore,  computationally  based  property  predictions  were 
obtained  to  support  consideration  of  2-azidocycloalkanamine  synthesis.  Included  are  DFT-based 
characterizations  of  trans  and  cis  ADMCPA  equilibrium  structures  and  their  gas-phase  heat  of 
formation.  A  limited  set  of  equilibrium  AMCBA  and  ACPA  conformers  is  also  characterized. 

In  addition,  an  approach  to  estimating  the  condensed-phase  densities  and  heats  of  vaporization 
for  azide-functionalized  compounds  based  on  molecular  dynamics  (MD)  simulations  with 
halogen  analogs  is  proposed  and  validated,  then  employed  to  predict  condensed-phase  properties 
for  the  ADMCPA  isomers.  The  property  predictions  are  subsequently  employed  in  ballistic 
performance  calculations  that  provide  a  comparison  of  ADMCPA  isomer  specific  impulses  (Isp) 
and  density  specific  impulses  (D*Isp)  relative  to  those  for  the  CINCH  fuels  and  MMH.  Finally, 
a  suggestion  for  synthesizing  czis- ADMCPA  is  proffered. 


2.  Computational  Quantum  Chemistry  Methods 


The  Gaussian98  (G98)  suite  of  quantum  chemistry  codes  [5]  was  employed  to  identify 
equilibrium  ADMCPA,  AMCBA,  and  ACPA  structures  and  characterize  their  normal  modes. 
The  level  of  theory  employed — ^DFT  utilizing  the  B3LYP  exchange-correlation  functionals  [6-8] 
and  6-311+-i-G(d,p)  atomic  orbital  basis  set  [9-11]  on  a  G98-defined  “ultrafine”  grid — was 
chosen  based  on  the  success  achieved  with  it  in  characterizing  DMAZ  [4].  In  addition  to  the 
structural  characterizations,  estimates  of  gas-phase  heats  of  formation  at  298  K  (AH/ (298))  were 
sought  for  use  in  predicting  ballistic  performance.  Unfortunately,  the  well-established  approach 
for  obtaining  AH/298)  estimates  via  the  G2  method  [12]  proved  computationally  intractable, 
necessitating  the  implementation  of  an  alternate  method.  The  results  presented  in  this  report  are 
derived  from  the  heats  of  reaction  (AHr(O))  calculated  for  “isodesmic”  reactions  at  0  K. 
(Isodesmic  reactions  are  those  in  which  the  total  number  of  each  type  of  bond  are  equal  in 
reactants  and  products;  the  success  of  the  method  assumed  to  be  associated  with  the  cancellation 
of  errors  achieved  through  equating  product  and  reactant  bonds  [13].)  In  an  attempt  to  identify  if 
the  method  generated  systematic  errors,  a  variety  of  reactions  involving  the  formation  of 
molecules  of  interest  and  xH2  (x  =  1  or  2)  from  reactants  with  accepted  AH/298)  values  were 
constructed  and  their  AH/298)  predictions  compared.  The  zero-point  corrected  energies  and 
enthalpy  differences  at  298  and  0  K  (H(298)-H(0))  for  the  reactants  and  reference  elements 
needed  to  derive  AH/298)  from  AH/O)  were  completed  at  the  same  level  of  theory  as  that 
employed  to  characterize  the  aliphatic  amine  azides  of  interest.  The  author  also  notes  that 
though  “accepted”  H(298)-H(0)  values  are  reported  in  the  Joint  Army-Navy-NASA-Air  Force 
(JANNAF)  thermochemical  tables  [14]  for  the  reference  elements  (Na,  Ha,  and  C(=solid))  and  a 
small  subset  of  the  reactants  (NH3  and  CH4),  for  internal  consistency,  the  H(298)-H(0)  values 
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determined  via  the  quantum  chemical  calculations  are  employed  in  all  cases  except  C(solid). 
(H(298)-H(0)  values  for  C(solid)  are  not  calculable  via  the  theoretical  methods  employed.) 


3.  Computational  Quantum  Chemistry  Results 


3.1  ADMCPA 

Figures  2  and  3  show,  respectively,  the  equilibrium  rra«i'-ADMCPA  and  m- ADMCPA 
structures  identified  through  this  study.  The  geometric  parameters  for  these  conformers  are 
provided  in  Table  1,  and  Table  2  presents  their  zero-point  corrected  energies,  dipole  moments, 
and  rotational  constants.  Beyond  their  classification  as  trans  and  cis  forms,  the  conformers  are 
distinguished  by  the  (dihedral  angle)  orientations  of  the  (CH3)2N-  and  -N3  groups  about 
(CH3)2N-(C3H4)  and  (C3H4)-N3  bonds,  and  the  labeling  scheme  was  chosen  such  that 
alphabetically  adjacent  conformers  are  nominally  transformed  into  one  another  by  a  dihedral 
angle  rotation  about  one  of  these  bonds.  The  scheme  also  groups  conformers  based  on  their 
H10-C9-N1 1-C12  angle’s  proximity  to  one  of  three  values.  For  rrcn^-ADMCPA,  these  angles 
are  -67°  (t-A  and  r-B),  61°  (t-C  and  r-D),  and  -165°  (t-E  and  r-F).  Conformers  t-A  and  r-B  are 
observed  to  be  more  than  4  kcal/mol  lower  in  energy  than  the  other  conformers,  indicating  that 
the  -67°  H10-C9-N1 1-C12  angle  is  energetically  favored.  (In  this  configuration,  the  amine  lone 
pair  points  toward  and  directly  over  the  -(C3H4)-  group.)  Two  different  N2-N3-C4-C6  dihedral 
angles  are  also  observed,  87°  (t-A,  t-D,  and  t-E)  and  - 165°  (r-B,  t-C,  and  r-F),  but  for  conformer 
pairs  where  this  is  the  only  primary  difference — (t-A,  r-B),  (t-C,  r-D),  and  (r-E,  t-F) — the  energy 
differences  are  small  (<0.2  kcal/mol).  The  six  trans  conformers  identified  are  all  of  the 
structures  that  can  be  constructed  fi'om  combinations  of  the  geometric  preferences  observed. 
Based  on  a  Boltzmann  distribution  calculated  fi'om  trans  confonner  zero-point  corrected 
energies,  it  is  expected  that  t-A  (58%)  and  t-B  (41%)  populations  will  dominate  a  gas-phase 
ADMCPA  sample  at  room  temperature. 

The  cw-ADMCPA  conformers  have  energetically  favored  H10-C9-N11-C12  and  N2-N3-C4-C6 
angles  similar  to  those  seen  in  the  rranj-ADMCPA  conformers,  but  the  correspondences  between 
conformer  energies  and  these  angles  are  slightly  different.  Perhaps  the  most  noteworthy 
difference  is  the  1.0  kcal/mol  difference  between  conformers  c-A  and  c-B.  This  difference  is 
presumably  related  to  the  fact  that  c-A  has  a  (stabilizing)  azide-amine  interaction  similar  to  that 
observed  in  the  lowest  energy  DMAZ  conformer  (Figure  1).  In  DMAZ,  the  conformer  with  this 
configuration  is  also  lower  in  energy  than  otherwise  similar  conformers  by  1.0— 1.3  kcal/mol. 
Therefore,  it  is  expected  that  c-A  will  prove  to  be  the  lowest  energy  cw- ADMCPA  conformer. 
The  six  cis  conformers  identified  are  all  of  the  structures  that  can  be  constructed  fi'om 
combinations  of  the  geometric  preferences  observed.  Based  on  a  Boltzmann  distribution 
calculated  from  cis  conformer  zero-point  corrected  energies,  it  is  expected  that  c-A  (82%)  and 
c-B  (18%)  populations  will  dominate  a  gas-phase  cw- ADMCPA  sample  at  room  temperature. 
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Figure  3.  Equilibrium  cis-ADMCPA  conformers.  (H8  of  conformers  c-B  and  c-C  is  hidden  in  the 
perspectives  shown.) 


Table  1 .  Geometric  parameters  for  equilibrium  ADMCPA  conformers  as  computed  via  B3LYP/6-3 1 1++G(d,  p). 
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1.091 

1.091  1.091 

1.093  1.093 

1.105  1.105 


172.9  172.9 


1.092  1.093  1.093 

_  _0„[  __  105  ^  ^  _ 

Simple  Angle  (in  Degrees) 


173.1  173.0  172.9 


115.8 

02.8 

20.2 


1.092  1.093 

1.106  1.101 


172.2 


116.7 

117.0  1 

1  116.6  1 

1  116.2 

1  116.5  1 

117.2  1 

i  117.2  1 

1  111.5  1 

1  111.5  1 

1  111.3 

llimiini 

116.2  111.3 


Table  2.  Relative  zero-point  corrected  energies,  dipole  moments,  and  rotational  constants  of 
equilibriiun  ADMCPA  conformers  as  con5)uted  via  B3LYP/6-31 1++G(d,p). 


Conformer  Designation 

AE 

(kcal/mol) 

u 

(Debye) 

Rotational  Constants 
(GHz) 

trans  conformers  \ 

t-A 

a 

2.56 

4.81 

0.76 

0.68 

t-B 

0.2 

2.53 

3.21 

0.93 

0.76 

t-C 

4.4 

2.53 

2.85 

0.99 

t-D 

4.4 

2.40 

4.07 

0.79 

0.73 

t-E 

4.8 

2.88 

5.03 

0.73 

0.71 

t-F 

5.0 

2.98 

3.31 

0.85 

0.79 

cis  conformers  \ 

c-A 

5 

2.95 

2.32 

1.31 

1.15 

c-B 

1.0 

2.87 

2.69 

1.14 

0.88 

c-C 

4.3 

1.94 

2.28 

1.28 

1.02 

c-D 

4.7 

1.74 

2.92 

0.99 

0.92 

c-E 

5.4 

2.87 

2.50 

1.14 

0.96 

c-F 

6.2 

2.64 

4.09 

0.86 

0.78 

“  Zero-point  corrected  energy  is  —415.401984  Hartrees. 
^  Zero-point  corrected  energy  is  -415.3991 13  Hartrees. 


The  relative  energies  of  the  lowest  energy  trans-  and  ci5- ADMCPA  conformers  also  deserve 
note.  That  is,  despite  c-A  having  the  stabilizing  amine-azide  interaction  found  in  DMAZ,  it  is 
still  1.8  kcal/mol  higher  in  energy  than  the  lowest  energy  fraws- ADMCPA  structure  (where  this 
interaction  is  precluded).  The  accuracy  of  the  energy  calculations  is  assumed  to  be  from  1  to 
3  kcal/mol,  so  die  true  energy  ordering  is  uncertain.  But  it  may  be  found  that  the  two  isomers 
yield  different  (overall)  reaction  exothermicities  in  an  otherwise  similar  system. 


With  19  atoms  and  Cs  symmetry,  all  ADMCPA  conformers  have  51  normal  modes.  Frequency 
values  and  (integrated)  infrared  (IR)  intensities  for  transitions  to  all  modes  of  each  conformer  are 
provided  in  Tables  3  and  4.  Lacking  experimental  mid-IR  absorption  spectra  with  which  to 
compare,  the  results  are  considered  to  have  limited  value  at  this  time,  and  so  are  not  discussed  in 
detail.  However,  they  do  indicate  that  if  ADMCPA  molecules  are  synthesized,  the  acquisition  of 
mid-IR  absorption  spectra  in  the  range  from  400  to  900  cm“^  will  reveal  which  structures  are 
present  in  the  sample.  Simulating  mid-IR  absorption  spectra  {I(v))  by  modeling  vibrational  band 
contours  (due  to  rotational  line  structure)  with  a  Lorentzian  lineshape. 


J 


y!  7t 

-yy+r'’ 


(1) 


where  Aj  is  the  integrated  absorption  for  a  transition  to  the  normal  mode  and  /(equal  10  cm"*) 
is  the  lineshape  width  at  half-maximum.  Figure  4  shows  that  the  presence  of  cw- ADMCPA  will 
be  identifiable  by  the  presence  of  a  feature  at  about  590  cm“*,  and  that  c-A  and  c-B  will  be 
identified/distinguished  by  the  concomitant  appearance  of  features  near  either  650  cm“^  or  700 
cm“'.  The  presence  of  a  feature  at  447  cm“',  on  the  other  hand,  will  identify  the  presence  of 
fra«5-ADMCPA. 


8 


Table  3.  Harmonic  vibrational  frequencies  and  integrated  IR  intensities  for  rra«j:-ADMCPA  conformers. 


Mode 

Conformer  | 

1  ^-A 

1 _ _ 

1  ^-c 

1 _ to _ 

1  /-E 

1 _ _ 

V 

A 

y 

A 

V 

1  ^ 

V 

A 

V 

rr- 

A 

1 

33 

0.0 

34 

0.2 

35 

31 

0.0 

33 

0.1 

2 

85 

0.4 

89 

0.5 

56 

56 

1  0.2 

49 

0.3 

3 

136 

1.1 

107 

0.4 

103 

1.3 

121 

0.6 

124 

108 

0.0 

4 

159 

0.8 

192 

1.2 

194 

0.1 

154 

3.1 

179 

0.6 

209 

2.3 

5 

225 

1.2 

222 

1.2 

227 

0.9 

e^h 

WKBM 

0.7 

225 

0.2 

6 

242 

0.3 

243 

249 

0.3 

IHSi 

0.4 

252 

0.3 

7 

278 

4.7 

270 

1.2 

266 

2.4 

279 

0.6 

1.5 

258 

1.4 

8 

308 

1.5 

328 

1.0 

327 

5.1 

327 

3.1 

1  292 

3.0 

323 

2.9 

9 

389 

3.6 

389 

10.6 

391 

8.0 

9.8 

381 

11.2 

10 

406 

4.9 

401 

5.6 

400 

1.8 

4.1 

426 

1.9 

11 

447 

2.1 

446 

1.7 

459 

2.1 

1.3 

464 

1.8 

458 

1.6 

490 

492 

1.9 

506 

5.5 

508 

2.2 

505 

5.9 

504 

4.3 

554 

1  9.1 

549 

8.4 

544 

5.8 

555 

8.4 

552 

7.4 

14 

654 

13.5 

642 

20.4 

645 

HEQI 

659 

15.6 

639 

■eh 

15 

■EBI 

785 

6.3 

762 

8.2 

756 

HBOi 

740 

4.9 

753 

IHEEI 

16 

1.3 

829 

7.1 

816 

6.7 

832 

7.9 

HSI 

818 

17 

36.4 

929 

34.5 

903 

HSB 

892 

18 

952 

9.9 

934 

931 

3.8 

932 

8.4 

915 

llllllllllll^ 

918 

8.4  1 

19 

969 

6.5 

977 

9.4 

961 

12.1 

942 

5.5 

966 

975 

20 

18.0 

1032 

15.3 

1027 

13.5 

1024 

1012 

19.6 

16.7 

21 

6.8 

1055 

9.8 

■ESH 

16.0 

1059 

14.5 

1040 

22.9 

1057 

22.2 

22 

MtiW 

16.7 

1079 

11.5 

■ESH 

5.8 

1070 

11.1 

1061 

15.3 

1060 

9.8 

23 

■ESI 

0.8 

1094 

10.2 

33.0 

1087 

11.9 

■eh 

ee^b 

5.4 

24 

nil 

2.3 

■IIEH 

7.5 

nil 

28.6 

■EEl 

edh 

23.2 

25 

MUM 

4.3 

■EH 

4.2 

■IBH 

5.8 

eeiei 

6.8 

1115 

5.2 

moQii 

26 

WIE-U 

1.4 

■EH 

2.4 

■EEQEI 

4.7 

■EE3I 

1.3 

1151 

5.9 

1155 

■■Bi 

27 

1186 

17.5 

1184 

16.0 

1174 

10.5 

1175 

8.8 

1175 

10.6 

H!^I 

H^^l 

28 

1223 

4.5 

5.2 

■B3I 

15.3 

■yjcj^E 

18.7 

20.8 

EBBH 

7.5  1 

29 

1243 

20.3 

■■gi 

10.8 

13.8 

■EH 

10.1 

hsh 

30 

1311 

6.9 

8.4 

13.3 

1314 

10.4 

EQEII 

34.8 

1310 

■^bb 

31 

1314 

38.4 

113.8 

1325 

195.5 

ebhh 

172.7 

1326 

MWaiM 

32 

1331 

210.2 

1331 

164.0 

1350 

40.6 

^ed 

■eh 

39.3 

1354 

51.9 

33 

1413 

4.1 

1406 

28.3 

1416 

■eei 

■E^H 

1.6 

1419 

13.2 

34 

1441 

1440 

0.0 

1442 

1442 

0.7 

0.4 

1441 

0.3 

msm 

4.9 

1465 

6.1 

1473 

1472 

2.3 

HH 

2.5 

EBH 

msm 

7.4 

1483 

7.1 

1484 

EH 

1484 

10.6 

EEHI 

21.1 

37 

9.7 

1485 

21.4 

■E^H 

28.2 

1486 

1493 

13.7 

1488 

19.4 

38 

6.2 

1497 

9.6 

■£^■1 

5.1 

1500 

7.1 

1499 

11.6 

1499 

8.1 

39 

■EH 

mm 

■^El 

25.0 

1502 

22.6 

15.6 

■eh 

HQBI 

40 

10.4 

msm 

msm 

■BQI 

16.6 

1518 

17.3 

1517 

15.2 

■eh 

hee^b 

41 

699.9 

2232 

626.8 

2228 

648.9 

2229 

686.9 

MHiM 

42 

49.3 

43.8 

64.4 

2911 

75.0 

EB^I 

43 

2929 

149.3 

2929 

EI^H 

^^E 

97.8 

Emi 

116.8 

2931 

120.4 

EH 

SljjeEl 

48.2 

■I^B 

51.2 

hsssh 

KlOi 

MM 

lEOi 

■eh 

33.6 

3056 

32.4 

BtWiM 

■ES 

3067 

26.0 

3067 

25.3 

38.6 

■BQEI 

IESH 

H3HI 

22.7 

EEOii 

21.2 

eh 

HQEl 

EEBI 

EHEi 

EE^B 

IHBI 

BBITilMH 

31.7 

48 

3105 

37.9 

3105 

31.1 

3124 

3.3 

3122 

2.9 

3127 

5.6 

3130 

5.8 

49 

3129 

8.9 

3132 

8.1 

3130 

22.7 

3128 

27.2 

3136 

4.9 

3135 

4.2 

50 

3134 

8.5 

3135 

8.6 

3136 

11.7 

3134 

9.2 

3141 

16.2 

3139 

19.3 

51 

3222 

4.2 

3225 

3.5 

3219 

5.1 

3218 

5.6 

3215 

6.6 

3218 

5.6 

Notes:  v=  cm“  y4  -  kilometers  per  mole. 
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Table  4.  Harmonic  vibrational  frequencies  and  iategrated  IR  intensities  for  cw-ADMCPA  conformers. 


Mode 

1  Conformer 

c-A 

r— 

1  c-C 

1  c-D 

1  c-E 

1  c-F 

V 

A 

V 

A 

A 

V 

A 

V 

A 

V 

A 

1 

50 

0.1 

32 

0.0 

0.1 

27 

0.1 

26 

0.1 

28 

0.3 

2 

85 

0.0 

77 

0.3 

69 

0.8 

75 

1.5 

IBEEB 

3 

1.1 

151 

3.9 

EIH! 

^BSi 

BUI 

beei 

4 

191 

1.0 

166 

1.4 

172 

1.4 

189 

5.9 

eh 

1.7 

213 

5 

227 

1.9 

227 

0.2 

214 

1.1 

213 

0.8 

1.1 

233 

6 

241 

240 

0.4 

228 

0.5 

225 

1.2 

ESSI 

257 

0.9 

7 

257 

■KfH 

259 

0.7 

282 

7.7 

280 

2.5 

^■Ql 

265 

1.5 

8 

Km 

5.6 

ESI 

5.1 

^El 

IEH 

3.8 

eeb 

1.2 

9 

1.3 

ESI 

5.0 

ESI 

eee 

BEE! 

2.2 

E^l 

BBS 

10 

409 

1.6 

HBH 

1  403 

eh 

3.2 

1.1 

1  420 

^L2J 

11 

496 

2.8 

■nQi 

kb 

3.2 

12 

537 

11.2 

545 

8.2 

548 

7.9 

553 

8.5 

1  551 

^BS 

13 

592 

3.8 

589 

5.4 

588 

21.9 

605 

574 

9.2 

BUB 

14 

701 

17.0 

656 

656 

6.7 

628 

662 

20.1 

bsb 

E^l 

15 

5.0 

753 

5.8 

718 

737 

7.5 

IIIIIIIIIIIIIIIII^QI^ 

16 

14.5 

808 

0.8 

1  782 

1  7.8 

808 

HE^I 

781 

4.4 

ieb 

17 

838 

12.6 

839 

IEH 

831 

7.4 

829 

3.7 

1.3 

18 

39.7 

IEH 

l^gl 

eei 

918 

933 

15.9 

931 

24.3 

19 

1.7 

IEH 

6.0 

1  965 

wmam 

955 

964 

6.0 

IEH 

7.9 

ebe 

eh 

3.6 

eh 

■bd 

1023 

6.5 

21 

HQH 

1  1058 

7.1 

1iT3M 

e^b 

eh 

5.0 

eh 

1039 

18.1 

22 

20.5 

TiVZ-B 

10.8 

MBlrM 

14.7 

1073 

1061 

28.8 

23 

6.5 

IBM 

eh 

19.0 

ITif 

9.3 

1092 

1093 

9.8 

24 

1117 

4.0 

1116 

4.2 

1116 

3.3 

3.8 

1115 

3.3 

5.4 

25 

1130 

1.0 

1129 

0.5 

1126 

in 

4.1 

1119 

1125 

1.4 

26 

1137 

2.5 

IEH 

4.9 

lltM 

1.7 

1159 

1.8 

HIEM 

1158 

14.8 

27 

16.1 

IEH 

14.4 

ww 

8.5 

1173 

5.7 

MIW 

1177 

12.6 

28 

1201 

1213 

11.3 

MUIM 

8.1 

1214 

23.0 

11.7 

29 

1242 

IEEH 

21.5 

19.4 

19.6 

13.6 

BkSE 

HBI 

30 

1314 

111111^3^^1 

1314 

5.2 

MItVX-M 

■MM 

41.0 

■MM 

31 

1361 

1337 

44.3 

EiSkl 

1350 

34.8 

■K 

45.1 

1338 

32 

WItVAM 

1368 

124.6 

eh 

111.3 

1378 

33 

1397 

24.7 

BEJiSJII 

eei 

besb 

■EIM 

9.2 

1425 

2.7 

34 

1441 

0.5 

1441 

0.2 

hese 

EwB 

0.6 

■EBM 

1.2 

2.3 

35 

3.5 

9.0 

eh 

2.6 

1472 

3.2 

eh 

8.5 

6.6 

36 

8.4 

5.7 

eh 

10.2 

1480 

13.9 

eh 

6.9 

1485 

10.7 

37 

1484 

16.9 

1485 

18.1 

1483 

9.9 

1489 

21.5 

1490 

19.5 

38 

1497 

6.1 

1497 

10.4 

BEkE 

8.0 

9.6 

BUthI 

7.0 

3.8 

39 

15.8 

1498 

11.3 

3.5 

1.9 

BLiB 

2.6 

■KiM 

8.7 

40 

1514 

1523 

17.9 

Mfcfef 

18.6 

1518 

■QB 

■MM 

15.8 

41 

EEI 

2231 

636.2 

2233 

753.1 

42 

Elfl 

lESII 

65.4 

2920 

iH 

43 

ElE 

lE^I 

■skliki 

77.8 

2928 

KKWH 

BEI 

44 

IESH 

HSDI 

3050 

E^l 

ESkI 

42.6 

3058 

29.0 

3056 

33.0 

45 

Ei^l 

El^l 

■aiSM 

memm 

28.5 

3068 

34.7 

46 

■aiM 

30.8 

ehi 

eh 

lE^I 

■eM:M 

30.1 

Bkkiil 

27.4 

3098 

29.9 

47 

eh 

21.8 

Eli^ll 

20.0 

EkB 

14.5 

EliB 

15.5 

eh 

12.2 

48 

eh 

18.1 

eim 

27.9 

eezi 

BtlBM 

eh 

10.8 

EK 

8.5 

49 

eh 

5.8 

3117 

18.4 

29.9 

3112 

33.1 

eh 

11.1 

eee 

3.4 

50 

eh 

9.2 

3130 

BlitM 

14.1 

3131 

11.9 

eh 

15.2 

■KM 

18.5 

51 

3210 

7.8 

3220 

_ isj 

3218  1 

5.9 

3218 

5.8 

3217  1 

6.9 

3213  1 

7.7 

Notes:  v  ~  cm'  ^  ;  A  =  kilometers  per  mole. 
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Figure  4.  Simulated  mid-IR  absorption  spectra  for  low  energy  rra«j-ADMCPA 
and  cis-ADMCPA  structures. 

3.2  AMCBA  and  ACPA 

The  equilibrium  AMCBA  and  ACPA  structures  identified  through  this  study  are  shown  in 
Figures  5  and  6,  respectively.  The  conformers  in  these  limited  sets  were  obtained  via  geometry 
optimizations  started  with  structures  expected  (based  on  the  ADMCPA  characterizations)  to  be 
energetically  favored.  Shielding  configurations  in  trans  isomers  were  also  sought.  Complete 
geometric  parameters  and  normal  mode  details  for  the  structures  in  these  figures  are  provided  in 
the  Appendix.  Briefly,  as  in  DMAZ  and  cw- ADMCPA,  the  cis  isomers  of  AMCBA  and  ACPA 
have  configurations  with  interacting  amine  and  azide  groups  that  are  1.0  kcal/mol  lower  in 
energy  than  otherwise  similar  conformers.  Thus,  a  population  dominated  by  a  shielded 
configuration  is  expected  for  the  cis  forms.  As  for  tranj- AMCBA  and  trans- A,  the 
increasing  flexibility  of  their  rings  allows  closer  proximity  of  the  azide  and  amine  groups  than 
in  fran^-ADMCPA;  the  N2-N(amine)  distance  measuring  their  separation  decreases  from 
4.29  A  in  fraw^-ADMCPA  [riB]  to  3.71  A  in  ?ran5-AMCBA  and  3.16  A  in  trans-AC? A.  The 
N2-N(amine)  separation  in  trans- AC? A  is  very  close  to  the  3.07-A  N2-N(amine)  separation  in 
the  lowest  energy  DMAZ  structure,  limiting  the  author’s  interest  in  further  characterization  of 
this  molecule. 
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Figure  6.  An  equilibriiun  ACPA  conformer.  (Arrow  points  to 
atom  hidden  in  the  perspective  shown.) 

3.3  Gas-Phase  Heats  of  Formation 

Table  5  shows  the  gas-phase  AHj(298)  estimates  obtained  for  various  alkanes,  aliphatic  amines, 
CINCH  compounds,  and  the  ADMCPA  isomers  based  on  calculated  AHr(0)  values  for  isodesmic 
reactions.  However,  one  of  the  concerns  with  such  estimates,  i.e.,  an  inability  to  identify  and 
correct  systematic  error  [13],  raises  a  question  about  their  reliability  and  requires  comment.  The 
audior  attempted  to  address  this  issue  by  constructing  isodesmic  reactions  in  a  consistent  manner, 
i.e.,  forming  the  molecule  of  interest  and  H2  from  reactants  with  accepted  AH/298)  values.  For 
alkanes  and  aliphatic  amines  with  accepted  AH/298)  values,  this  method  yields  reasonably 
accurate,  but  systematically  high  (0-3.5  kcal/mol),  estimates.  This  suggests  that  the  method’s 
values  can,  in  fact,  be  rationalized  and  was  the  basis  for  extending  the  method  to  obtain  AH/298) 
estimates  for  the  azide-functionalized  molecules  in  the  table. 

Another  general  concern  in  predicting  AH/298)  via  theoretical  results  for  isodesmic  reaction 
schemes — namely,  the  necessity  of  good  experimental  heat  of  formation  values  for  all  (assumed) 
knowns — proves  to  be  an  issue  for  the  azide-functionalized  fuels  and  also  requires  comment. 
First,  the  bonds  in  the  azido  group  are  unique  in  nature — ^i.e.,  they  are  not  reliably  represented  as 
single,  double,  or  triple  bonds,  or  by  other  “common”  bonds — ^and  thus  the  isodesmic  schemes 
must  include  an  azido  compoimd  as  a  reactant.  Unfortunately,  the  only  simple  azide  compound 
widi  a  reliable  heat  of  formation  value  is  hydrozoic  acid  (HN3).  Heats  of  formation  for  methyl 
azide  and  ethyl  azide  would  have  been  valuable,  but  published  experimental  values  could  not  be 
found,  and  there  are  indications  that  values  obtained  with  the  G2  method  are  from  2  to 
3  kcal/mol  too  low  [15].  Thus  the  ability  to  compare  values  from  different  isodesmic  schemes 
and  identify  systematic  errors  is  limited.  Based  on  the  author’s  “best  estimates”  for  the 
AH/298)' s  of  methylazide  and  ediylazide  and  an  experimentally  obtained  AH/298)  for 
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Table  5.  Gas-phase  AHf(298)  estimates. 


Compound 

ZPE 

H(298)-H(0) 

AHKO) 

AHK298) 

Diff. 

Isodesmic  Reaction  Scheme 

{a.u.) 

JANNAF 

(kcal/mol) 

Calc 

(kcal/mol) 

Exp* 

(kcal/mol) 

Calc 

(kcal/mol) 

Exp** 

(kcal/mol) 

Calc 

(kcal/mol) 

C(solid) 

0.25 

0.0 

0.00 

H2 

-1.169506 

2.02 

2.07 

0.0 

0.00 

N2 

-109.554124 

2.07 

2.07 

0.0 

0.00 

Ammonia 

-56.548462 

2.40 

2.40 

-9.3 

Methane 

-40.489427 

2.39 

-16.0 

-17.90 

Cyclopropane 

-117.850067 

2.71 

17.0 

12.75 

Hydrazoic  acid 

-164.814684 

2.62 

71,8 

70.30 

Ethane 

-79.782286 

BOB 

-18.9 

1.3 

CH4+CH4-^Ethane+H2 

Butane 

-158.374256 

4.23 

BliUI 

-26.7 

3.5 

2  Ethane^Butane+Hz 

Methylamine 

-95.830097 

2.75 

-1.8 

-1.3 

-5.0 

0.5 

CH4+NH3->MeNH2+H2 

Ethylamine 

-135.129229 

3.37 

-5.6 

-4.7 

CH4+MeNH2^EtNH2+H2 

-5.5 

Ethane+NH3->>EtNH2+H2 

Cyclopropylamine 

-173.202101 

24.1 

24.6 

19.0 

(CH2)3+NH3->CpNH2+H2 

Dimethylamine 

-135.117679 

2.5 

-4.50 

-2.9 

CH4+MeNH2^Me2NH+H2 

-213.715623 

4.95 

-8.8 

-5.3 

-13.8 

3.5 

EthanfrbEtNH2->Et2NH+H2 

-174.408125 

4.02 

1.4 

3.4 

^0291 

-3.7 

2.0 

CH4+Me2NH->Me3N+H2 

Methyl  azide 

-204.103331 

3.40 

75.5 

75.5 

(73) 

72.4 

(-0.6) 

CH4+HN3-^MeN3+H2 

Ethyl  azide 

-243.402659 

mi^^9 

71.1 

71.1 

66.4 

(-1.6) 

Ethane+HN3->EtN3+H2 

1-azidobutane 

-321.995173 

55.9 

ButanefHN3->BuN3+H2 

MMAZ 

5.56 

86.5 

78.8 

'  CH4+EtNH2+HN3->>MMAZ+2H2 

DMAZ 

lom 

78.8 

Me2NH+Ethane<-HN3-^DMAZ+2H2 

CPAZ 

PIHEliRMd 

(CH2)3+EtNH2+HN3-^CPAZ+2H2 

/ra«5-ADMCPA 

-415.401984 

miBsm 

106.7 

Me2NH+(CH2)3+HN3-^ADMCPA+2H2 

cw-ADMCPA 

-415.399113 

6.56 

117.5 

108.4 

Me2NH+(CH2)3+HN3^ADMCPA+2H2 

'  “Experimental”  AH|(0)  values  are  the  sum  of  experimental  (exp)  AHr(298)  values  and  calculated  (calc)  H(298)-H(0)  values. 

’’  Values  in  parentheses  are  estimates  proposed  by  the  author  based  on  consideration  of  results  published  by  Rogers  and  McCafferty  [15]  and  Lee  et  al.  [16]. 


1-azidobutane  [16],  the  aliphatic  azide  AH/298)  values  appear  to  be  undeipredicted.  But  since 
these  undeipredictions  are  coupled  with  (expected)  slight  overpredictions  for  the  (isodesmic) 
formation  of  secondary  and  tertiary  amines  from  primary  and  secondary  amines,  the  author 
expects  that  the  Table  5  AH/298)  estimates  for  the  aliphatic  amine  azides  from  these  schemes 
will  prove  to  be  reasonably  accurate.  An  experimental  value  for  DMAZ  has  been  established 
[2,  3]  and  indicates  that  the  theoretically  based  prediction  is  approximately  3  kcal/mol  (5%)  too 
high.  However,  given  the  lack  of  corroborating  evidence  for  the  experimental  value,  this  issue  is 
considered  an  open  one. 


4.  Molecular  Mechanics  Methods  and  Results 


In  addition  to  the  gas-phase  properties  characterized,  the  author  was  interested  in  estimating 
condensed-phase  densities  and  heats  of  vaporization  for  each  ADMCPA  isomer.  Coupled  with 
the  gas-phase  heats  of  formation  determined,  such  property  predictions  can  be  employed  in 
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calculations  to  predict  a  fuel’s  performance  in  a  propulsion  system.  Having  been  employed 
successfully  by  Bunte  and  Sun  for  predicting/estimating  the  density  of  nitrate  esters  [17],  an 
approach  for  obtaining  condensed-phase  properties  based  on  the  Molecular  Simulations,  Inc. 
(MSI,  now  Accelrys,  Inc.)  suite  of  molecular  modeling  codes  and  its  COMPASS  force  field  [18] 
was  desired.  However,  COMPASS  force  field  parameters  for  azide  atom  types  had  not  been 
developed  at  the  time  of  the  study,  so  the  use  of  “azide  analogs”  in  the  simulation  routines  had  to 
be  considered.  The  properties  of  notional  azide  compounds  are  commonly  predicted  based  on 
their  halogen  analogs,  and  a  specific  approach  for  predicting  densities  was  suggested  by  the  work 
of  Lee  et  al.  [16],  who  found  that  the  (experimentally  determined)  densities  of  6>-azidoalkanes 
are  nearly  identical  to  those  of  their  monochloride  analogs.  In  addition,  these  researchers  found 
that  the  heats  of  vaporization  of  cu-azidoalkanes  are  close  to  those  of  their  monobromide  analogs. 

To  establish  that  the  correspondences  observed  by  Lee  et  al.  would  be  found  computationally, 
MD  simulations  with  the  monochloride  and  monobromide  analogs  of  DMAZ,  CPAZ,  MMAZ, 
and  1-azidobutane  were  set  up  and  run.  The  simulations  were  performed  on  cubic  three- 
dimensional  periodic  cells  containing  1080-1100  atoms  (56-72  molecules)  constructed  by  the 
Amorphous  Cell  module  of  Insight  II.  (Amorphous  Cell  builds  cells  according  to  a  modified 
Markov  process,  incorporating  rotational  isomeric  state  probabilities  modified  to  account  for 
nonlocal  [long-range]  interactions  [19,  20].)  Eight  different  cells  with  an  initial  target  density  of 
0.9  g/cm^  were  built  for  each  compound,  the  length  of  cube  edges  being  ~25  A  for  monochloride 
compounds  and  ~23  A  for  monobromide  compounds.  A  group-based  cut-off  method  with  tail 
correction  was  employed  to  evaluate  nonbond  interactions.  The  cut-off  method  assumes  that  the 
radial  distribution  fimctions  converge  to  rmity  beyond  the  cut-off  distance.  The  cut-off  distance 
was  specified  to  be  10  A  for  all  compounds  except  the  CPAZ  analogs,  where  a  10.5-A  cutoff  was 
employed. 

The  simulations  consisted  of  three  stages,  the  dynamics  in  all  cases  being  modeled  with  Verlet 
velocity  integration  [21]  and  Andersen  temperature  control  [22].  Berendsen  pressure  control 
[23]  was  employed  for  constant  pressure  simulations.  The  first  stage  was  a  30,000-step, 

1-fs/step,  constant  volume  and  temperahire  (NVT)  “pre-equilibration”  simulation  run  to  relieve 
large  stresses  (inadvertently)  introduced  by  the  cell  packing  procedure.  The  second  stage  of  the 
procedure  was  a  30,000-step,  1-fs/step  constant  pressure  and  temperature  (NPT)  simulation  that 
allowed  the  cell  to  equilibrate.  The  final  stage  was  a  50,000-step,  1-fs/step  NPT  simulation 
during  which  data  was  recorded. 

Table  6  compares  the  densities  of  the  monochloride  analogs  of  DMAZ,  CPAZ,  MMAZ,  and 
1-azidobutane  vs.  experimentally  determined  values  for  the  corresponding  azides.  The 
computationally  determined  density  for  1-chlorobutane  is  observed  to  be  nearly  identical  to  its 
experimentally  determined  value  and  that  of  1-azidobutane.  (Given  that  the  COMPASS  force 
field  parameters  for  the  atom  types  relevant  to  1-chlorobutane  were  probably  developed  with  a 
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Table  6.  Condensed-phase  property  estimates. 


Compound 

Density 

AHv(298) 

AHK298) 

MD 

Exp 

Exp 

Calc 

Exp 

Calc 

Exp 

(g/cm^ 

(g/cm^) 

(kcal/mol) 

(kcal/mol) 

(kcal/mol) 

(kcal/mol) 

(cal/g) 

(cal/g) 

1-azidobutane® 

9.17 

46.5 

50.7 

470 

510 

l-chlorobutane* 

0.882  ±0.003 

0.886 

7.8 

8.01 

1-bromobutane* 

7.5 

8.79 

2-azido-N-Tnethylethanamine  (MMAZ)** 

0.970 

69 

690 

2-chloro-N-methylethanamine 

0.980  ±0.003 

9.6 

2-bromo-N-methylethanamine 

9.2 

2-azido-N,N  dimethylethanamine  (DMAZ)** 

0.933 

9.4 

69 

66 

580 

2-chloro-N,N  dimethylethanamine 

0.955  ±0.004 

9.5 

2-bromo-N,N  dimethylethanamine 

8.8 

2-azido-N-cyclopropylethanamine  (CPAZ)** 

0.993 

94 

750 

2-chloro-N-cyclopropylethanamine 

0.981  ±0.006 

9.1 

2-bromo-N-cyclopropylethanamine 

9.7 

■■ 

fro«5-2-chloro-N,N-dimethylcyclopro- 
panamine  (/ra/i.s-ADMCPA  analog) 

■ 

J  enantiomer  only 

0.946  ±0.006 

1^1 

Racemic  mixture 

9.5 

97 

cw-2-azido-N,N-dimethylcyclopro- 
panamine  (cw-ADMCPA  analog) 

■ 

1  enantiomer  only 

0.992  ±0.006 

Racemic  mixture 

1.000  ±0.005 

10.6 

98 

790 

Monomethylhydrazine  (MMH)® 

0.883 

9.6 

13.1 

284 

•  Experimental  data  are  from  Lee  et  al.  [16]. 

**  Experimental  data  are  from  Thompson  et  al.  [2,  3]. 

®  Experimental  data  are  as  reported  in  the  Schmidt  [1]  compendium. 


training  set  consisting  of  <2>-alkylchlorides,  the  good  agreement  was  expected.)  The 
computationally  determined  densities  for  the  monochloride  analogs  of  MMAZ,  DMAZ,  and 
CPAZ  are  also  observed  to  be  in  reasonable  agreement  with  experimentally  determined  values 
for  the  corresponding  azide.  The  approach  predicts  that  the  density  of  ^raMJ-ADMCPA  will  be 
0.945  g/cm^,  a  value  that  is  slightly  lower  than  that  predicted  for  DMAZ  (0.955  g/cm^).  The 
density  of  cw-ADMCPA  is  predicted  to  be  1.00  g/cm^,  a  value  that  is  even  higher  than  that 
predicted  for  CPAZ  (0.98  g/cm^). 

Table  6  also  presents  estimates  for  the  heats  of  vaporization  of  the  monochloride  and 
monobromide  analogs  of  DMAZ,  CPAZ,  MMAZ,  and  1-azidobutane.  It  is  observed  that  the 
values  for  the  monochloro  and  monobromo  analogs  of  DMAZ  and  1-azidobutane  are  similar  to, 
but  lower  than,  the  values  observed  experimentally  for  the  corresponding  azides.  Since  the  heats 
of  vaporization  determined  for  the  monochloride  and  monobromide  analogs  are  similar,  the 
effort  needed  to  perform  a  separate  set  of  simulations  with  bromide  analogs  will  not  be 
imdertaken  in  the  future.  Being  a  byproduct  of  simulations  necessary  to  predict  densities,  the 
heats  of  vaporization  derived  from  monochloride  analogs  are  “free,”  and  thus  provide  an  easily 
obtained  reference  point.  In  the  present  case,  however,  it  is  observed  that  most  of  the  values  in 
the  table  cluster  near  DMAZ’s  experimentally  determined  value  of  9.4  kcal/mol,  and  this  value 
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was  employed  by  the  author  in  calculating  all  of  the  condensed-phase  heats  of  formation  given  in 
Table  6.  The  potential  error  associated  with  this  approach  is  expected  to  be  less  than 
±2  kcal/mol. 


5.  Ballistic  Properties 


With  estimates  for  condensed-phase  heats  of  formation  and  densities,  it  is  possible  to  calculate 
and  compare  the  specific  impulses  (Isp)  and  density  impulses  (D*Isp)  that  MMAZ,  DMAZ, 
CPAZ,  and  the  ADMCPA  isomers  will  produce  when  combined  with  an  oxidizer.  Table  7 
compares  optimized  Isp  values  from  NASA-Lewis  rocket  calculations  [24]  for  a  system  in  which 
the  oxidizer  is  an  IRPNA  “mixture.”  (The  oxidizer  mixture  is  one  employed  by  AMCOM  for 
such  comparisons.)  It  is  observed  that  the  optimized  Isp  values  for  the  aliphatic  amine  azides  are 
nearly  identical  to  each  other  and  less  than  1%  lower  than  the  value  obtained  for  MMH. 
Moreover,  with  higher  densities  than  MMH,  their  D*Isp’s  are  predicted  to  be  higher. 


6.  Synthesis  Possibilities 


The  author  has  no  (relevant)  experience  in  the  art  of  synthetic  organic  chemistry,  but  wanted  to 
at  least  attempt  to  address  the  possibility  of  synthesizing  ADMCPA.  A  literature  search  of 
synthesis  routes  for  aliphatic  amine  azides  revealed  a  large  number  of  possibilities,  with  a 
relatively  recent  review  of  the  preparation  of  azide  compounds  written  by  Schriven  and  Tumball 
being  particularly  instructive  [25].  Of  the  approaches  identified  via  this  survey,  the  most 
intriguing  to  the  author  was  suggested  by  the  work  of  Benalil  et  al.  [26].  Based  on  synthesizing 
desired  1,2  aminoazides  by  treating  “appropriate”  1,2  aminoalcohols  with  triphenylphosphorous 
(PPhs),  diethyl  azodicarboxylate  (DEAD),  and  HN3 — i.e.,  the  (modified)  Mitsunobu  reaction 
[27] — ^this  approach  should  be  readily  testable  because  an  appropriate  1,2  aminoalcohol  for 
ADMCPA,  namely  rra«5-2-hydroxy-N,N-dimethylcyclopropanamine,  has  already  been 
synthesized  by  Armstrong  and  Cannon  [28].  Since  the  modified  Mitsunobu  reaction  converts 
alcohols  to  azides  through  inversion  of  configuration,  successful  conversion  of  fra«5'-2-hydroxy- 
N,N-dimethylcyclopropanamine  via  this  approach  will  probably  produce  cz5-ADMCPA. 


7.  Discussion 


The  superior  performance  of  hydrazine  and  its  methylated  derivatives  compared  to  other 
hypergols  can  generally  be  attributed  to  three  characteristics;  short  ignition  delays,  high 
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Table  7.  NASA-Lewis  thermochemical  code  Isp  and  D*Isp  predictions  for 
various  fuels.® 


Compound 

O/F 

I 

SP 

%A 

Throat 

(s) 

Exit 

(s) 

D*Isp‘’ 

(Ibfs/in^ 

2.0 

112.9 

284.8 

2.2 

111.4 

MMAZ 

2.4 

109.9 

284.4 

+10 

2.6 

281.0 

2.8 

107.1 

276.7 

3.0 

272.0 

2.0 

281.1 

2.2 

112.9 

283.3 

DMAZ 

2.4 

284.4 

9.8 

+7 

2.6 

284.4 

2.8 

■nsMjB 

282.8 

3.0 

107.8 

279.5 

2.0 

113.9 

281.3 

2.2 

113.2 

283.6 

CPAZ 

2.4 

112.0 

284.8 

2.6 

110.8 

284.9 

10.1 

+10 

2.8 

283.9 

3.0 

281.1 

2.0 

2.2 

113.4 

284.6 

cw-ADMCPA 

2.4 

112.3 

285.6 

2.6 

285.7 

10.3 

+12 

2.8 

284.5 

3.0 

108.5 

281.7 

2.0 

282.1 

2.2 

■1^ 

284.2 

2.4 

msam 

MSSBM 

9.8 

+7 

2.8 

3.0 

281.5 

2.0 

115.6 

284.7 

115.0 

286.8 

MMH 

2.4 

288.0 

2.6 

HESHI 

288.1 

9.2 

— 

110.9 

285.7 

280.8 

®  Results  for  “infinite-area”  combustor  problem  with  2000-psia  chamber  pressure  and 
equilibrium  chemistry.  The  subsonic  area  ratio  is  15.0,  and  the  supersonic  area  ratio  is 
14.175. 

^  For  consistency,  in  all  cases  (except  MMH),  the  density  predictions  obtained  via  the 
molecular  modeling  computations  are  employed  to  calculate  D*Isp. 
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densities,  and  high  specific  impulses,  the  shortcoming  of  most  alternatives  being  long  ignition 
delays.  Given  that  aliphatic  azides  are  not  hypergolic  and  that  the  ignition  delays  of  aliphatic 
amines  are  much  longer  than  those  of  hydrazine  and  its  derivatives,  the  short  ignition  delays 
produced  by  (the  aliphatic  amine  azide)  DMAZ  are  surprising,  and  identifying  the  mechanisms 
that  imderlie  this  phenomenon  could  be  useful  in  designing  fuels  with  even  shorter  ignition 
delays. 

Since  aliphatic  azides  are  not  hypergolic,  but  do  undergo  (explosive)  exothermic  unimolecular 
decomposition  if  raised  to  a  sufficiently  high  temperature,  and  aliphatic  amines  are  hypergolic 
but  “slow,”  in  the  absence  of  detailed  knowledge  about  DMAZ’s  molecular  structure,  a  “working 
hypothesis”  was  suggested.  That  hypothesis  was  that  exothermic  proton  transfer  from  nitric  acid 
to  the  amine,  which  is  more  nucleophilic  than  the  azide,  was  jump-starting  the  azide 
decomposition,  and  the  azide  decomposition  was  accelerating  the  otherwise  slow  secondary 
reactions  involved  in  pure  amine  oxidation.  This  hypothesis  was  consistent  with  the  observation 
of  ammonium  nitrate  salts  in  combustion  residue  [3],  and  observed  differences  in  the 
performance  2-azidoethanamines  with  different  amine  substituents,  shorter  ignition  delays  found 
for  fuels  whose  amines  would  be  expected  to  be  stronger  bases.  Not  surprisingly,  this  hypothesis 
led  AMCOM  to  focus  on  synthesizing  fuels  with  more  a  more  nucleophilic  amine. 

The  genesis  of  this  study  was  the  finding  that  in  DMAZ’s  lowest  energy  configuration,  the  azide 
group  sterically  shields  the  lone  pair  electrons  associated  with  the  amine  nitrogen  [4].  Given  the 
experimentally  indicated  importance  of  proton  transfer  to  this  site,  the  author  initially  assiuned 
that  shielding  would  reduce  a  compound’s  reactivity  and  therefore  contribute  to  increasing 
ignition  delays.  Thus,  his  interest  was  in  identifying  compounds  that  would  prevent  shielding. 
Through  the  course  of  this  study,  however,  several  observations  led  the  author  to  consider  the 
possibility  that  the  amine-azide  interaction  in  the  shielding  configuration  is  actually  the  enabling 
mechanism.  First,  given  the  fact  that  close  to  50%  of  DMAZ’s  amine  sites  will  be  shielded  from 
proton  attack  at  room  temperature  and  that  DMAZ  is  competitive  with  Aerozine  50,  it  would  be 
imprudent  to  ignore  this  possibility.  Second,  though  (aliphatic)  azides  are  not  known  as 
hypergols,  the  azide  group  has  two  lone  pair  sites  and  nitric  acid  hydrogen  bonds 
(exothermically)  to  these  sites  in  the  gas  phase.  Calculations  for  an  as  yet  incomplete  study  to 
characterize  DMAZ-nitric  acid  complexes  indicate  that  hydrogen-bonding  to  azide  sites  is  only 
slightly  less  than  the  exothermic  than  the  hydrogen-bonding  at  amine  sites.  In  addition,  though 
enhancing  the  amine’s  basicity  would  almost  certainly  reduce  ignition  delays,  there  is  the 
potential  that  any  increase  in  amine  basicity  will  be  offset  by  an  enhancement  of  the  amine-azide 
interaction  and  an  increase  in  shielding.  (Such  an  effect  is  suggested  by  preliminary  results  from 
ongoing  theoretical  characterizations  of  CPAZ  and  MMAZ.)  If  enhancement  of  the  azide 
group’s  basicity  through  the  azide-amine  interaction  is  the  enabling  mechanism,  then  cis- 
ADMCPA  could  yield  very  short  ignition  delays.  The  author  hopes  that  the  promise  of  the 
ADMCPA  isomers  that  is  indicated  in  the  computational  characterizations  presented  is  sufficient 
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to  pursue  their  synthesis,  as  experimental  testing  with  them  should  clarify  the  mfluence  of  the 
azide-amine  interaction  and  the  azide’s  basicity  on  the  ignition  process. 

As  a  final  note,  the  possibility  of  using  a  CH2  group  to  bridge  the  ethyl  group  carbons  of  DMAZ 
was  investigated  because  of  the  obvious  smulanties  between  the  2-azidocycloalkanamines  and 
DMAZ  and  CPAZ,  but  the  bridge  can  presumably  be  made  with  groups  other  than  CH2,  NH  and 
O  being  two  possibilities.  A  stable  structure  based  on  the  former— an  aziridine— is  shown  in 
Figure  7.  However,  the  author  does  not  have  enough  knowledge  of  small  heterocyclic  organic 
compounds  to  even  guess  as  to  whether  their  properties  will  be  similar  to  DMAZ  or  CPAZ,  so  he 
made  no  attempt  to  further  characterize  such  alternatives.  Certainly,  property  estimates  for  such 
alternatives  could  be  calculated  if  considered  warranted. 


Figure  7.  An  equilibrium  3-azido-2-dimethylaminoaziridine  conformer. 


8.  Summary 


Fuels  based  on  2-azidocycloalkanamines  are  proposed  as  means  of  obtaining  hydrazine- 
alternative  hypergols  with  ignition  delays  shorter  than  those  observed  for  DMAZ.  The  geometric 
parameters  and  normal  modes  for  a  number  of  2-azidocycloalkanamines  with  C5HioN4 
stoichiometry  were  characterized  through  the  use  of  DFT,  and  MD  simulations  performed  with 
the  halogen  analogs  of  ADMCPA,  DMAZ,  CPAZ,  and  MMAZ  were  employed  to  predict  liquid 
densities  and  heats  of  vaporization.  Isp  and  D*Isp  values  predicted  based  on  the  DFT  and  MD 
results,  coupled  with  the  possibility  of  shorter  ignition  delays,  suggest  that  the  ADMCPA 
isomers  merit  further  consideration  as  alternatives  to  DMAZ  or  MMH  in  hypergolic  propulsion 
system  applications. 
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Appendix.  Details  From  2-Azido-N-methylcyclobutanamine  (AMCBA)  and 
2-Azidocyclopentanamine  (ACPA)  Theoretical  Characterizations 


This  Appendix  provides  tables  summarizing  the  theoretieal  characterizations  of  cis  and  trans 
forms  of  2-azido-N-methyl-cyclobutanamine  (AMCBA)  and  a  trans  isomer  of  2-azido- 
cyclopentanamine  (ACPA).  Table  A-1  presents  the  geometric  parameters  for  the  AMCBA 
conformers  shown  in  Figure  5,  and  Table  A-2  presents  their  relative  zero-point  corrected 
energies,  dipole  moments,  and  rotational  constants.  Table  A-3  presents  the  normal  mode 
frequencies  and  integrated  infrared  (IR)  intensities  for  transitions  to  these  modes  for  the 
AMCBA  structures.  Table  A-4  presents  the  geometric  parameters  for  the  ACPA  conformer 
shown  in  Figure  6,  and  Table  A-5  presents  its  normal  mode  frequencies  and  the  integrated  IR 
intensities  for  transitions  to  these  modes.  The  dipole  moment  for  the  ACPA  conformer  is 
3.94  Debye,  and  its  rotational  constants  are  3.07, 1.15,  and  0.88  GHz.  Its  zero-point  corrected 
energy  is  -415.448209  Hartrees. 
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Table  A-1.  Geometric  parameters  for  equilibrium  AMCBA  conformers  as  computed  via 
B3LYP/6-311-H-G(d,p). 


Conformers* 

Parameter 

c-A  1  c-B 

r-A  1  t-D 

1  Bond  Length  (in  Angstroms 

_ _ 

N1-N2 

1.136 

1.135 

1.135 

1.136 

1.228 

1.230 

1.230 

1.229 

1.475 

1.471 

1.469 

1.474 

C4-C6 

1.546 

1.546 

1.542 

1.553 

H5-C4 

1.095 

1.095 

1.097 

1.095 

C6-C9 

1.553 

1.553 

1.553 

1.553 

H7-C6 

1.089 

1.089 

1.090 

1,090 

H8«C6 

1.091 

1.091 

1.092 

1.092 

C9-C12 

1.550 

1.546 

1.570 

1.549 

H10-C9 

1.092 

1.092 

1.094 

1.092 

H1I-C9 

1.091 

1,091 

1.090 

C12-C4 

1.567 

1.577 

1.549 

1.546 

H13-C12 

1.105 

1.103 

1.094 

1.103 

N14-C12 

1.447 

1.441 

1.442 

1.449 

H15-C14 

1.016 

1.012 

1.015 

1.015 

C16-C14 

1.464 

1.459 

1.462 

1.463 

H17-C16 

1.093 

1.094 

1.094 

1.092 

H18-C16 

1.092 

1.092 

1.092 

1.092 

H19-C16 

1.103 

1.103 

1.102 

1.102 

1  Simple  Angles  (in  Degrees 

_ 

N3-N2-N1 

173.0 

172.4 

173.1 

173.7 

C4-N3-N2 

116.8 

116.8 

116.8 

116.5 

H5-C4-N3 

109.9 

110.2 

110.4 

110.7 

C6-C4-N3 

108.8 

108.9 

115.7 

00 

bo 

H7-C6-C4 

115.3 

115.5 

117.3 

117.3 

H8-C6-C4 

111.5 

111.6 

110.3 

110.9 

C9-C6-C4 

89.1 

89.1 

88.5 

87.7 

H10-C9-C6 

111.3 

111.0 

112.0 

111.9 

H11-C9-C6 

117.9 

118.3 

117.6 

117.6 

C12-C9-C6 

88.1 

88.6 

87.8 

88.6 

H13-C12-C9 

108.7 

109.6 

109.3 

110.1 

N14-C12-C9 

118.4 

118.9 

117.2 

117.5 

H15-N14-C12 

110.0 

110.9 

110.4 

109.4 

C16-N14-C12 

113.6 

115.8 

115.0 

114.2 

H17-C16-N14 

109.6 

109.9 

109.2 

109.6 

H18-C16-N14 

109.3 

109.5 

109.4 

109.4 

H19-C16-N14 

113.9 

113.7 

114.2 

113.7 

1  Dihedral  Angles  (in  Degrees)  1 

C4-N3-N2-N1 

168.9 

-166.3 

174.4 

- 179.3 

H5-C4-N3-N2 

-56.1 

58.5 

-36.8 

42.4 

C6-C4-N3-N2 

173.7 

-171.1 

164.2 

-88.7 

H7-C6-C4-N3 

24.2 

-24.7 

96.5 

101.5 

H8-C6-C4-N3 

149.7 

-150.4 

-29.2 

-25.0 

C9-C6-C4-N3 

-96.1 

95.7 

-142.0 

-137.8 

H10-C9-C6-C4 

92.1 

-92.5 

-92.1 

-92.5 

H11-C9-C6-C4 

-139.9 

139.7 

140.6 

139.6 

C12-C9-C6-C4 

-19.4 

19.2 

19.4 

H13-C12-C9-C6 

-88.5 

88.7 

89.0 

89.2 

N14-C12-C9-C6 

140.9 

-141.1 

-141.5 

-139.6 

H15-N14-C12-C9 

-172.5 

-30.2 

49.2 

-50.9 

C16-N14-CI2-C9 

64.5 

-158.7 

76.4 

175.0 

H17-C16-N14-C12 

-64.7 

-65.2 

-62.7 

-64.4 

H18-C16-N14-C12 

177.6 

177.4 

179.7 

177.6 

H19-C16-N14-C12 

56.2 

55.8 

58.2 

56,0 

‘  Labels  refer  to  Figure  5. 
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Table  A-2.  Relative  zero-point  corrected  energies,  dipole  moments,  and  rotational  constants  for 
equilibrium  AMCBA  conformers  as  computed  via  B3LYP/6-31 1++G(d,p). 


Conformer  Designation^ 

AE 

(kcal/mol) 

(Debye) 

Rotational  Constants 
(GHz) 

trans  conformers 

t-A 

_ b 

3.54 

2.24 

1.15 

0.80 

t-B 

0.2 

2.41 

1.10 

0.80 

CIS  conformers 

c-A 

1.2 

3.25 

2.10 

1.43 

0.92 

c-B 

_ 

3.38 

2.03 

1.57 

0.99 

^  See  Figure  5. 

^  Zero-point  corrected  energy  is  -415.414753  Hartrees. 
^  Zero-point  corrected  energy  is  -415.410703  Hartrees. 
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Table  A-3.  Harmonic  vibrational  frequencies  and  integrated  IR  intensities  for  equilibrium  AMCBA 
conformers. 


Mode 

1  Conformer*  | 

1  c-A 

1  c-B 

1  /-A 

_ _ 

V 

1  ^ 

V 

1  ^ 

V 

A 

V 

L_^l _ 

1 

58 

64 

35 

0.2 

27 

2 

97 

1.3 

79 

1.0 

84 

77 

1.9 

3 

105 

0.5 

104 

1.0 

106 

116 

1.7 

4 

173 

1.3 

175 

0.9 

172 

0.8 

147 

5 

192 

0.4 

183 

0.6 

197 

0.4 

6 

225 

0.3 

251 

1.4 

236 

1.0 

235 

7 

298 

1.7 

273 

2.9 

278 

4.2 

332 

1.9 

8 

367 

6.3 

3.2 

368 

355 

8.5 

9 

439 

5.7 

417 

4.0 

449 

3.1 

419 

5.5 

10 

519 

4.6 

530 

13.8 

491 

3.2 

474 

2.7 

11 

559 

16.0 

548 

5.1 

8.8 

12 

647 

HBSHi 

633 

552 

4.8 

17,6 

13 

699 

2.0 

667 

21.1 

680 

39.8 

7.9 

14 

751 

8.3 

695 

10.5 

735 

738 

84.3 

15 

777 

59.0 

755 

1.0 

818 

19.5 

811 

16 

840 

4.3 

837 

7.9 

822 

3.3 

4.9 

17 

905 

S.l 

3.1 

5.5 

924 

3.1 

18 

937 

1.3 

933 

1.8 

933 

1.3 

927 

1.3 

19 

950 

8.4 

962 

2.8 

991 

12.0 

997 

38.2 

20 

1004 

6.4 

996 

4.3 

14.5 

1007 

0,1 

21 

1041 

18.1 

■ESH 

9.6 

16.3 

1049 

2.9 

22 

1076 

9.4 

8.6 

1066 

7.5 

1078 

11.9 

23 

1116 

5.6 

1114 

1.4 

■OH 

12.4 

■OH 

1.2 

24 

00 

1.9 

1143 

4.0 

mssm 

■DEH 

25 

1173 

13.3 

1169 

24.6 

1162 

38.2 

26 

1201 

22.5 

1187 

11.4 

1199 

19.1 

1.8 

27 

1212 

1.1 

1235 

3.5 

1211 

2.9 

1212 

2.6 

28 

1246 

0.8 

1248 

0.7 

1236 

0.6 

1242 

0.2 

29 

1266 

1.1 

1262 

2.2 

1252 

0.3 

1255 

1.8 

30 

1268 

5.1 

1272 

9.8 

1265 

1263 

2.0 

31 

1294 

6.4 

1286 

9.5 

1280 

15.3 

1284 

2.0 

32 

138.4 

1332 

123.6 

1329 

114.4 

1329 

135.3 

33 

HESSHl 

65.9 

1342 

58.1 

1356 

64.1 

1353 

77.0 

34 

34.4 

1383 

22.3 

1394 

65.7 

1402 

55.0 

1459 

1463 

0.1 

1458 

1.9 

1460 

3.0 

36 

1475 

7.3 

4.3 

1483 

3.8 

1474 

3.1 

37 

1480 

2.6 

1476 

8.3 

1491 

14.6 

1482 

38 

1493 

8.8 

1493 

8.4 

1496 

6.3 

1493 

39 

1496 

6.4 

1495 

6.5 

9.1 

1502 

10.1  1 

40 

1518 

11.2 

1517 

12.7 

1522 

27.1 

1518 

BOH 

41 

2231 

609.7 

583.4 

2233 

684.5 

2230 

42 

2907 

71.8 

36.2 

2953 

mmm 

41.4 

43 

2938 

114.2 

2944 

3009 

mmm 

103.0 

44 

47.7 

3026 

45.3 

18.2 

3041 

14.6 

45 

3059 

7.7 

3052 

25.2 

23.3 

3056 

34.8 

46 

3062 

3053 

34.4 

3062 

35.5 

47 

3068 

31.2 

32.2 

25.3 

3062 

18.2 

48 

23.1 

24.7 

3097 

25.6 

3102 

29.7 

49 

3111 

16.3 

24.8 

3102 

26.2 

3105 

50 

3127 

27.0 

24.7 

3120 

34.7 

51 

3506 

0.2 

3556  1 

3.3 

3516 

HHEHI 

3521  1 

^hk^hI 

‘  See  Figure  5. 

Notes:  y  -  cm"*;  j4  =  kilometers  per  mole. 
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Table  A“4.  Geometric  parameters  for  an  equilibrium  ACPA 

conformer  as  computed  via  B3LYP/6-31 1++G(d,p). 


Parameter* 


Bond  Length  (in  Angstroms) 


H8-Ce 

1.093 

C9-C6 

1.557 

HIO-C 

1.092 

Hll-C 

1.092 

C12-C 

H13-C 

C15-C 

] 

1.543 

C15-C 

] 

1.542 

H16-C 

1 

1.096 

N17-C15 


H18-N17 

H19-N17 


N3-N2-N1 


C4-N3-N2 

H5-C4-N3 


C6-C4-N3 


H7-C6-C4 


H11^9-C6 

C12-C9-C6 


Simple  Angles  (in  Degrees) 


172.2 


H8-C6-C4 

108.2 

C9-C6-C4 

105.6 

H10-C9-C6 

111.8 

H14-C12-C9 

113.2 

C15-C12-C9 

105.1 

H16-C15-C4 

112.0 

N17-C15-C4 

111.9 

H18-N17-C15 

111.6 

H19-N17-C15 

112.8 

Dihedral  Angles  (in  Degrees) 

C4-N3-N2-N1 

168.3 

H5-C4-N3-N2 

-38.2 

C6-C4-N3-N2 

-160.0 

H7-C6-.C4-N3 

84.4 

H8-C6-C4-N3 

-33.3 

C9-C6-C4-N3 

-152.2 

H10-C9-C6-C4 

-114.7 

HI  1-C9-C6-C4 

126.9 

C12-C9-C6-C4 

6.7 

H13-C12-C9-C6 

-98.4 

H14-CI2-C9-C6 

141.7 

C15-C12-C9-C6 

18.7 

HI6-C15-C4-N3 

46.4 

N17-C15-C4-N3 

-70.7 

H18-N17-C15-C4 
|H19~N17-C15-C4 
*  See  Figure  6. 
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Table  A-5.  Harmonic  vibrational  frequencies  and  integrated  IR 
intensities  for  an  ACPA  conformer  as  computed  via 
B3LYP/6-31H-+G(d,p)." 


Mode 

y 

A 

1 

40 

0.2 

2 

61 

0.5 

3 

114 

0.3 

4 

180 

1.6 

5 

231 

29.0 

6 

274 

1.1 

7 

354 

7.4 

8 

387 

11.4 

9 

469 

4.0 

10 

533 

3.5 

11 

556 

9.3 

12 

603 

1.5 

13 

691 

11.8 

14 

769 

11-5 

15 

796 

75.4 

16 

815 

43.1 

17 

910 

4.4 

18 

936 

16.7 

19 

956 

3.6 

20 

988 

7.0 

21 

1028 

6.5 

22 

1072 

9.5 

23 

1092 

0.5 

24 

1120 

6.4 

25 

1136 

10.6 

26 

1195 

1.4 

27 

1244 

1.9 

28 

1265 

0.5 

29 

1283 

2.0 

30 

1317 

0.5 

31 

1319 

0.6 

32 

1328 

2.9 

33  i 

1335 

124.9 

34 

1344 

2.5 

35 

1376 

65.3 

36 

1428 

12.0 

37 

1490 

0.8 

38 

1497 

5.2 

39 

1519 

7.1 

40 

1660 

34.4 

41 

2236 

655.3 

42 

2980 

20.6 

43 

3012 

19.5 

44 

3034 

31.1 

45 

3048 

26.7 

46 

3058 

35.5 

47 

3070 

24.6 

48 

3083 

17.3 

49 

3102 

51.0 

50 

3499 

0.5 

51 

3586 

3.2 

‘  See  Figure  6. 

Notes:  V  =  cu'  *;  ^4  =  kilometers  per  mole. 
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